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We derive an alternative representation for the relativistic non-local kinetic energy operator and 
we apply it to solve the relativistic Salpeter equation using the variational sine collocation method. 
Our representation is analytical and does not depend on an expansion in terms of local operators. 
We have used the relativistic harmonic oscillator problem to test our formula and we have found that 
arbitrarily precise results are obtained, simply increasing the number of grid points. More difficult 
problems have also been considered, observing in all cases the convergence of the numerical results. 
Using these results we have also derived a new representation for the quantum mechanical Green's 
function and for the corresponding path integral. We have tested this representation for a free 

"j^" ■ particle in a box, recovering the exact result after taking the proper limits, and we have also found 

that the application of the Feynman-Kac formula to our Green's function yields the correct ground 
, state energy. Our path integral representation allows to treat hamiltonians containing non-local 

£\j ■ operators and it could provide to the community a new tool to deal with such class of problems. 



o 
o 



PACS numbers: 03.65.Go,02.70. Jn,11.15.Tk 



i— » 

m 

I. INTRODUCTION 

m 
> 

■ The appearance of non-local operators in the relativistic extensions of the Schrodingcr equation poses a serious 
| challenge both to analytical and numerical calculations. However, the inclusion of relativistic effects is crucial for 
. example in the study of meson phenomenology, where the Bcthe-Salpetcr equation (BSE) provides the correct theo- 
retical tool to describe relativistic bound states. Replacing the kernel in the BSE with an instantaneous local potential 
one obtains a relativistic Schrodingcr equation, which is also known as "spinlcss Salpeter equation" (SSE). In such 
case the hamiltonian operator is typically given by 



H = vV + m 2 + V(r) . (1) 



Oh, From a technical point of view, the inclusion in the Hamiltonian of the relativistic kinetic energy operator, 
\/p 2 + m 2 , complicates the solution of the problem because of its non-local nature. The great phenomenological 
relevance of the SSE has motivated in the past twenty years many efforts to solve this equation, either using ana- 
lytical or numerical techniques. Early work on this subject is contained for example in [l[ 0, H, 0, H[ . The method 
k>( | described in which allows one to obtain a matrix representation of the non-local kinetic energy operator, has also 
j_J ■ been used recently Q in conjunction with the Lagrange mesh method. In this case, the method is particularly appeal- 
ing since it does not require the evaluation of the matrix elements of the potential, but rather only the specification of 
thepotential on the grid points. This feature is also shared by the so called sine collocation methods (see for example 
0, [![), which could also be used straightforwardly together with the method of [J]: however we do not consider this 
possibility since we are rather interested in developing a completely new approach, as it will soon become clear. 

In a series of recent works, Lucha and collaborators have been able to obtain precise upper and lower bounds to the 
eigenvalues of the RSE, 0, [13, [ll|, H2]: these bounds provide useful analytic or semi-analytic expressions and have 
been applied to a number of test problems (see for example 0). Another approach has been followed in 0, where 
an effective hamiltonian of non-relativistic form, which includes relativistic effects by means of parameters depending 
on the momentum, was constructed. 

Finally, the SSE was also studied in the case of the relativistic harmonic oscillator (RHO) in three recent papers, 
0,111, EH- In [3 the authors were able to obtain recurrence relations for the coefficients of the scries defining the 
cigenfunctions. The reader can also find useful the detailed bibliographic information contained in that paper. 

The main purpose of the present paper is to develop a completely new approach to the solution of the SSE: we 
will derive an analytical representation of the non-local kinetic energy operator and use it with the Variational Sine 
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collocation Method (VSCM) Q to obtain arbitrarily precise numerical results. We will use the RHO of for 
which semi-analytical results are available, to test our approach. The new representation that we have found for the 
relativistic kinetic energy operator has also been generalized to the calculation of the quantum mechanical Green's 
function and has allowed us to obtain a new formula, which differs from the standard formula of path integrals. 

The paper is organized as follows: in Section [IT] we describe the SSE for the RHO and obtain precise numerical 
solutions working in momentum space configuration (these results are then compared with the results of [HI); in 
Section imi we obtain an explicit analytical expression for the matrix elements of the non-local kinetic energy operator 
in terms of the so called "little sine functions" (LSF) recently studied in Q (these functions were first introduced by 
Schwartz in [I?) and later used by Baye in [Hj]); we use this representation in the VSCM working in coordinate space 
and reproduce the results previously obtained in momentum space; in Section IIVI we extend our method to obtain 
explicit analytical representation for the Green's function which holds for general potentials; finally in Section IVl we 
draw our conclusions and set the direction for future work. 



II. THE RELATIVISTIC HARMONIC OSCILLATOR 



The RHO problem corresponds to the case in which a spherical harmonic oscillator potential, V(r) 
in the Hamiltonian of eq. (TTJ). After a simple rescaling of the mass and of the energy the Schrodingcr equation can be 
cast in the form 



VP+^ + r 2 \ |V) =e |V) , (2) 

depending on a single parameter \i 1 . The advantage of considering this potential lies in the fact that the corresponding 
Schrodingcr equation in momentum space representation is local and it can therefore be attacked with standard 
techniques. In this case we are left with the equation 



A P + vV + V? *(P) = e *(P) , (3) 



where A„ = 



+ K 

The authors of [16( focus their analysis on the I = solutions and obtain a recurrence relation for the coefficients 
of the reduced radial wave function y(p) = y/Anp \P{p). The same problem has also been discussed by M. Trott in 
[l9| , using a basis of harmonic oscillator wave functions in momentum space and numerically computing in this basis 
the matrix element (i|y — d% + m^lj). 

We will use the VSCM of @, Q to obtain a numerical solution of the problem and then compare our solution 
with the results of [l6|. The VSCM uses sine functions (SF), defined on the real line, or "little sine functions" 
(LSF)H, a particular generalization of sine functions defined on finite intervals, to solve the Schrodinger equation by 
a collocation technique. Since the details of this method are clearly explained in 0, Q , we will here avoid mentioning 
all the technical details focusing on giving a more qualitative picture. 

SF and LSF are functions which are strongly peaked around a certain value and they fastly decay and oscillate, when 
moving away from this value. Under certain conditions they can be chosen to be orthogonal: such set of orthogonal 
functions is obtained by performing a replica of one function at the points where this function vanishes, thus effectively 
introducing a grid. Unless otherwise specified, cither by the Physics of the problem or by convention, the spacing 
of the grid is arbitrary and, if not carefully chosen, it strongly affects the precision of the results. In a recent paper 
Amore et al. have used an arbitrary scale factor as a variational parameter in the solution of the Schrodinger equation 
using a basis of simple harmonic oscillator wave functions [2^] : in that case it was shown that the optimal scale factor 
could be chosen to minimize the trace of the Hamiltonian matrix in the Hilbert subspace spanned by the N elements 
of the basis. The same principle, which was inspired by the Principle of Minimal Sensitivity (PMS) [2l| . was then 
used in 0, [f| using a sine collocation technique. As mentioned in the Introduction the sine collocation has the great 
advantage of avoiding the evaluation of matrix elements of the potential, which is simply "collocated" at the grid 
points; on the other hand the evaluation of the matrix elements of the non relativistic kinetic operator is also quite 
simple, since it involves matrices whose elements are obtained by collocating the derivative of a given sine function 
at the grid points. In this way one easily obtains a matrix representation for the Hamiltonian. 



1 We adopt the same notation used in 
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H ki = 4i >+S k iV(kh) . (4) 



For example in one dimension one has 

2m' 

where is the matrix obtained from the second derivative of the sine functions, h is the grid spacing. In the 
following we will set % = 1, so that h will not be confused with the Planck constant. 

The diagonalization of H allows one to obtain the eigenvalues (energies) and the eigenvectors (wave functions) of 
the problem (the number of these eigenvalues and eigenvectors being equal to the number of sine functions used) . 

Before trying to deal with the RHO in coordinate space, where it is non-local, we wish to use the VSCM to find 
numerical solutions in momentum space, where it is local. In this case machinery of VSCM applies straightforwardly 
and we are able to verify our claim of accuracy of our method. To allow a comparison with [lq we have used fi = 30 
and calculated the values of e — fi for /i = 30 using grids with different N. The optimal region in momentum space 
has been obtained by applying the PMS condition to the problem. Our results show that the method converges quite 
rapidly; as a matter of fact, for N = 50 we obtain e - \i = 0.386266042572445193517188444455 which has all the 
digits correct and agrees with the result shown in Table 1 of [l6j], although the latter contains only 6 digits. 

However we have also performed a more accurate test using the recurrence relations for the coefficients appearing 
in the series for the wave functions given in [l6| . We have extracted these coefficients by expanding around p = the 
wave function obtained with our method and we have compared them with the results obtained using directly the 
recurrence relation of eq. (8) of [l6[ ■ In this way we have obtained an indipendent confirmation of the accuracy of our 
results. 

The importance of the RHO lies in the unique possibility of having a complete control on the solutions of the 
problem, which can be calculated to any desired accuracy, despite retaining the non-local nature of the kinetic energy 
operator. In more general problems, when the potential is not quadratic in the coordinates, one cannot recover a 
local Schrodinger equation by working in the momentum space configuration. In these cases one possibility is to 
resort to a non-relativistic expansion in powers of p/ fi, which to lowest order provides the rest mass and the usual 
non-relativistic kinetic energy: 

p2 

v/p 2 + fi 2 » n + |- + . . . (5) 

However the hamiltonian operator obtained in this way contains derivatives of higher order. In such cases it is still 
possible to apply the VSCM in coordinate space to solve these equations, by using the matrices corresponding to the 
higher order derivatives (which can be computed quite easily). Although the implementation of this procedure poses 
no problem, we wish to show that in certain cases it can provide unexpected results, such as results which converge 
to wrong solutions as the higher order relativistic corrections are added. 

We use once again the RHO as our "guinea pig" and work in momentum space configuration using a potential 
obtained by expanding the kinetic term up to a give order in p 2 / fi 2 . The expression in eq. ([5]), for example, would 
represent the potential obtained working to order p 2 / n and corresponds to the usual non-relativistic harmonic 
oscillator. However, in the inclusion of the higher order terms we have to take into account that terms of order 
(p 2 /fi 2 ) 2n , with n integer, are negative, and therefore one always needs to work to order (p 2 /fj, 2 ) 2n+1 to have a 
spectrum bounded from below. We call V2 n +i(j>) the potential in momentum space corresponding to the expansion 
to order (p 2 / fj, 2 ) 2n+1 . For example to order (p 2 //z 2 ) 3 we have 

? 4 fi 

fr I \ i V P . P fa , 

Vsip) = fi+ 2;-8^ + i6^- (6) 

Once we substitute the non-local operator with its local expansion to a given odd order we can solve the corre- 
sponding local Schrodinger equation and thus obtain the energies and the wave functions. We have calculated the 
values of e — fj, for the ground state of the potential using N — 100 and two different values of fj, (fx — 30 and \i = 5). 
In the case corresponding to /i = 30 we have observed that the series converges to the exact result within 28 digit 
precision for V31 , while in the case corresponding to \i = 5 we have observed that the series does not converge to the 
exact result, providing its best approximation for Vg. 

We can give a simple explanation of this behavior: for p > (i the series in eq. ([5]) - where p is now a real number - 
diverges and therefore the potential obtained using this series, Voa(p) = lim, woo V2 n +i(p), corresponds to the original 
potential with an infinite wall located at p = [i. In the case \x = 30 the wave function is extremely small at p = \i 
and therefore the non-relativistic expansion is capable of providing highly (but not arbitrarily) accurate results; on 
the other hand, in the case fi = 5 the wave function is sizeable at p = \x and the non-relativistic expansion provides 
a very poor approximation. The situation would clearly get worse as /z is taken to be smaller. 

It should be stressed that in both cases the non-relativistic expansion docs converge, although not to the exact 
result, but rather to a value which depends on the scale \i. I am not aware if this is a well-known fact in the literature. 
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III. MATRIX ELEMENTS OF y/p 2 + p?- 

In this Section we want to describe a new way of calculating the matrix elements of the relativistic kinetic energy 
operator, which avoids the problems of the local non-relativistic expansion described in the previous Section. 

As mentioned in the Introduction, there are approaches in the literature to calculate the matrix elements of this 
operator, although they are numerical. For example, the method described in [J consists of 4 steps, i.e. the compu- 
tation of (K 2 )ij = (i\p 2 + /i 2 |j), the diagonalization of K 2 , the computation of the diagonal square root matrix and 
finally the computation of the square root matrix K. Because of the calculation of K is numerical this procedure 
would not be profitable in a variational scheme, where the PMS condition is always analytical. 

However, we will now show that it is possible to obtain once and for all an analytical representation of ^/p 2 + fi 2 and 
indeed of any non-local operator which is function of the momentum operator. Because our results will be analytical 
it will be possible to extend the VSCM to include relativistic non-local terms, in what we will call "Relativistic 
variational collocation method" (RVSCM). 

Let us now present our results. We will work in the following with the Little Sine Functions (LSF) of IS| , although 
similar results can also be obtained with the usual sine functions 2 . The LSF have been obtained in [8j using the 
orthonormal basis of the wave functions of a particle in a box with infinite walls located at x = ±L: 



= sin {^-jr{x + L) 



(7) 



A LSF is simply obtained as 



2L 

s k (h,N,x) = —^2ipn(x)ip n (y k ) 



N 



N 



(8) 



where y k = 2 j^ ± = kh. The LSF should be regarded as an approximate representation of a Dirac delta function. 
To simplify the notation we have introduced the grid spacing h = 2L/N. k is an integer which takes the values 
k = — N/2 + 1, —N/2 + 2, . . . , N/2 — 1, each corresponding to a different grid point. 

As shown in Q an explicit expression for the LSF can be obtained for even values of N in the form 



<(h N T )= 1 / sin((l + ^) l(x-kh)) cos((l + ^) f(x + kh)) 
After simple algebra one can also obtain the alternative expression 



(9) 



1 +N 

Sk (h,N,x)=-- Y, si 
n=— JV 



mr knir \ n^x 

sin [ 1 e 2L 

2 N 



(10) 



which is suited to calculate the action of ^Jp 2 + y 2 over a LSF: 



n=-N \ ' 



(11) 



Upon collocation on the grid, i.e. after setting x — > Xj, we have the matrix element 



K kj = V¥T7> s k (h, N, Xj ) = -I g (i)^ sin + *g) JQ) 



n=-N 



[i 2 e 1 n 



(12) 



2 As a matter of fact it is shown in § that for N — » oo and keeping h fixed the LSFs reduce to the SFs. 
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Notice that, despite its appearance, the expression above is real, as it would be evident expressing it in terms of 
trigonometric functions; however, we prefer to leave it in terms of plane waves. 

To check that eq. (|T2| is indeed the correct matrix element of the non-local kinetic energy operator we can calculate 
the matrix product (K 2 )ki = Ylj KkjKji 



N/2-1 +N +N 



(K\ t = y. k « k » = w E E E (»t 



+n 2 +2 



71l7T 



j'=-JV/2+l ni=-JV n 2 = -N 



~2~ + N 



sin 



n 2 7r jn 2 n 
2 N 



n 2 ir 



+ fi 2 e l " 



(13) 



2L J n V V 2L 

and compare it with the matrix (K 2 )ki = {p 2 + H 2 )ki, calculated either using the very same eq. (p~2|) with 



(ir) 2 + M 2 

We can define 



§f )~ + M 2 ~~ * (if)" + ^ or usm S the matrix for the second derivative to represent p (see 



N/2-1 

E,m . / n 2 7r jn 2 n 
e N sm ' 



j=-N/2+l 



N 



which for n 2 = ±ni takes the values 



C =-C = - — (-i) nl+1 



(14) 



(15) 



When C„i„2 is used inside eq. (fT3"|) it can be seen that only the terms corresponding to n 2 = ±ni contribute, whereas 
the remaining terms cancel out. In this case we obtain 



{K 2 



1 1:1 



1 

7p 



+N +N 



E E ( i ) ni+ri2 



m=-N n 2 = -N 



N 



Hit: kniir 
~ + N 



~2L 



n 2 Tr 
~2L 



/i 2 e l N 



„ ( i) + ) (^)ii,n 2 $m,— n 2 ) 



1 +N 



n=-N 



nix kmr\ 



\2l) + > r 



(16) 



This completes our proof, since we have precisely obtained the expression that we would have reached if we had used 
eq. (|12p directly with the squared operator. We have also verified numerically this result and confirmed its validity. 
As done for the non-relativistic Schrodinger equation we can obtain the Hamiltonian matrix 



H k i = K H + S k iV(kh) 



(17) 



For a fixed N the matrix elements Hu depend upon the arbitrary scale L: as discussed in 0, it is possible 
to obtain an optimal value for L by applying the principle of minimal sensitivity (PMS) to the trace of the N x N 
hamiltonian matrix TJv(L) = Tr [Hj\rxJv] : 



0. 



(18) 



Physically one can justify this condition by observing that the trace of the hamiltonian is an invariant under unitary 
transformations and therefore independent of L. On the other hand the trace of the truncated matrix does depend 
on L: one can therefore apply the PMS condition Q18p to minimize such dependence and thus obtain the scale where 
the problem is less sensitive to changes in L (normally this equation provides a single solution, although in the case of 
multiple solutions one would choose the solution corresponding to a flatter curve). Since (|18|) is an algebraic equation, 
the computational cost needed in solving it is quite limited. The scale L obtained in this way is then used inside 
eq. (fTT)) and the hamiltonian matrix is diagonalized, thus providing the N lowest eigenvalues and eigenfunctions. The 
error corresponding to choosing L according to the PMS is observed to be nearly optimal, leading to results which 
converge more rapidly as N increases. 

We have used eq. (jTTJ) in the case of the RHO with fi = 30 and /i = 5, i.e. the case which was previously studied in 
the momentum representation. The numerical results that we have obtained in this nonlocal representation display 
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a rate of convergence similar to the one observed in the local representation. In Fig. [T] we have plotted the difference 
|gW — e exac t\ as a function of the scale L and compared the results with the predictions obtained using the PMS 
condition eq. (fT5|) . Here is the energy for the ground state obtained using a given value of TV, whereas e exact is 
approximated with the energy obtained using N = 100. 

In the case /i = 5 where the non-relativistic expansion provides a quite poor approximation with just a four digit 
precision, the result obtained in the non-local representation has 12 correct digits for N = 50: [e— /x) = 0.91531941208. 

Although the results obtained studying the RHO are sufficient in our opinion to prove the efficiency and simplicity 
of the method that we are proposing in this paper, we wish to consider few examples which can better illustrate the 
power of our method. 

As a first example we have chosen the potential 



V[x) = Va 2 + . 



and the corresponding SSE (with x £ (0, oo)): 



\/p 2 + /i 2 + \J a 2 + x 2 ip(x) = e tp(x) 



(19) 



(20) 



Notice that for a ^ 1 and /i> 1 one recovers the standard harmonic oscillator (SHO). We have applied our method 
using a = [X = 1 and we have obtained (e — fi) = 2.075870921, using N = 100, where all the displayed digits are 
correct. 

As a second example that we wish to consider the SSE (with x £ (0, oo)): 



\J ffi + p 2 + /i 2 + \J a 2 + x 2 ip{x) = e ij){x) 



(21) 



Also in this case we have assumed \i = • 
the displayed digits are correct. 

As a final example we have considered the equation 



1 and we have obtained (e — /i) = 2.724199599, using N — 100, where all 



H log 



1 



P 

2fi 2 



+ V(x) \ ip(x) = E ip(x) 



(22) 



with V(x) = x 2 /2 and x £ (— oo, +oo). Notice that the equation obtained expanding to leading order in powers of 
p/fi is the standard Schrodinger equation for the harmonic oscillator. Despite the rather ugly form of the "kinetic" 
operator in the present case, the application of our method is straightforward and an expression for the matrix 
elements of the Hamiltonian is easily obtained: 



n=-N 



8fi 2 L 2 



— +5 kj V{x k ) 



(23) 
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We have solved this problem assuming jj, — 10 and we have obtained the ground state energy e = 20.04909725 using 
N = 30 (all the displayed digits are correct). 

Although we did not have any physical model in mind when we considered the last two examples, we are aware that 
modifications of the standard Schrodinger equation (not necessarily of this kind) appear in several areas of research: for 
example, the introduction of a minimal length uncertainty relation naturally leads to modified Schrodinger equations 
(see for example (22j|). 



IV. GREEN'S FUNCTIONS 



In this Section we want to show that the method that we have developed in the previous Section can also be applied 
to the calculation of matrix elements of other non-local operators. We will see that SF and LSF are a powerful 
tool, and actually they have already been used earlier in a different context in applications to Quantum Field Theory 
(QFT) (see for example [13). 

To start we consider an operator O = f(p): in such case, the results that we have obtained in the previous Section 
apply straightforwardly and one obtains 



n=-N x ' 

One example of application of this formula is the calculation of the Green's function for a free particle in a box, 
|x| < L, which is given by (see for example 0, HH) 

G^(y,0;x,t) = (x\e^\y) . 

The notation \x) indicates a state localized at a point x: in our formalism this is simply represented by a SF or a 
LSF with a peak at this point. Notice however that SF and LSF are not normalized to one, see 0, which means 
that we have an extra factor 1/y/h = \J N/2L for each SF or LSF, where h is the spacing of the grid. 

Therefore we can write: 

n=-N ^ ' 

Xk and Xj being points on the grid. In the limit of an infinitely dense grid one can switch from the discrete indices 
k,j to the continuum indices x,y: 

1 + N 2 2 

G ( %,0;M) = ~ £ ^e^-W + f ) (26) 

n=-N 

Notice that in the case of a free particle the exact result can be calculated with the standard path integral and 
reads 



G<o) (0j0;M)= /_£_ e *^ . (27) 



In our formalism we have 

+N 



G (o) (0,0;:M) = -j? £ i« +1 e" 1F " 1 ^ i sin (^) (28) 

-N 



1 + N / 2 2 2 2 2 

^ n=-N/2 



For N 3> 1 we can approximate this sum with an integral and thus obtain which reads 

' l , i 



1 / fi iM *? 

e 24 



2 V 2nit 



Erf f (j + |)((A r + l)^-2L jg) Erf f + «N-l)nt + 2L»x) 



30) 
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In the limit N — > oo we can approximate each error function to one and therefore obtain the exact expression (1271) . 
To further test this formula we can also use the Feynman-Kac formula to extract the ground state energy 

E = - lim - log / dx G(x, -it, x, 0) . (31) 

r — >oc T J 

In our case the integral appearing in (|3ip becomes a sum over all the grid points and therefore reads 



f 1 N/2-1 +N 
k=- N/2+1 n=-N 



So = -_Um _zlog<m >: >: sin(^ + ^) eA^Y] jH? y (32) 

Because we are taking the limit t — > oo, the exponential factor in this expression will be quite small unless ?i = ±1 
(the term n = vanishes). For this reason we can approximate 



- lim - log \ — V 

fc=-iV/2+l 



N/2-1 



i kir 
■ cos I — 



= -^M e ^} = w (33) 

which is indeed the exact ground state of the particle in the box. 

We can now consider the more general case in which the Hamiltonian contains a potential. In this case the formula 
of Section HVl cannot be applied directly because p and x do not commute. However we can use the Trotter product 
formula to write 0, [25[ 

/ .2 \N 

G(y,0;x,t) = lim (x\ I e~ l ^ t e~ i N' t ) \ y ) . (34) 

Af-»oo \ / 

We can use the completeness of the coordinate states, \x), to write 

G(y, 0; x, t) = \im^ £ (x\ Lr^e-^ \x ri )(x ri \ (V^V 4 ^ 1^) . . . (x rf , \ {e^^ e -^\ 1^35) 

where the indices span all the lattice. Notice the factors t/N, which correspond to having a "time-slicing", as usual 
in path integration; at each intermediate time each point on the grid can be reached. 
We define 

V^t) = h (x^e-^lx,) = e-^*% , (36) 
which is diagonal in the coordinates. We can now write the compact expression 

G(y,0;x,t) = lim V G (0) (x ri , t - e; x, t) V^(x ri , t) G^(x r . 2 , t - 2e; x ri , t - e) 

r 1 ,r 2 ,...,rf l 

V N {x r21 t) ... Vjf{x rft ,t) G<®(y,0]x rtt ,e) , (37) 
where e = t/N. This expression should be compared with the standard path integral expression [2^, [25| 

G(x,t;y,0) = lim f dx x . . . dx N ^ (JL.)"' 2 [^C 21 ^)'-^)] (38) 

where e = t/N. 

Our equation (|37|) can be also compared with the exact expressions for the TV-fold time sliced spacetime propagators 
which Crandall has obtained in eqs.(2.9) and (2.12) of [26[ using the standard representation. 

Our result of eq. (|37| provides a new representation of the path integral for quantum mechanical problems with a 
clear physical interpretation: the propagation of a particle sitting at y at time t = and reaching x at time t occurs 
moving at each time interval, e, from a point of the grid to another one in all possible ways (remember the sums over 
the grid points). G^°'(xi,t]Xj,t + e) represents the probability of going from a point Xi to another point Xj on the 
grid in a time interval e; at this point an interaction take place, through the the potential term V. We stress that 
our representation corresponds to a different way of discretizing space which also allows to deal with Hamiltonians 
containing non local operators (the SSE is one example) and could be an useful tool for problems which cannot be 
easily treated with the conventional formalism. 
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V. CONCLUSIONS 



In this paper we have derived an analytical expression for the non-local relativistic kinetic energy operator which 
appears in the Salpeter equation. This representation is exact in the limit of an arbitrarily fine grid and we have used 
it to solve the Salpeter equation for the Relativistic Harmonic Oscillator, where semi-analytical results arc available. 
We have found that our representation can be used together with the Variational Sine Collocation Method (VSCM) 
to provide arbitrarily precise results, with a strong rate of convergence. The most important result of this paper is the 
new representation for the quantum mechanical Green's function, which requires the evaluation of matrix elements of 
non-local operators. We have provided a general formula and we have explicitly tested it in the case of a free particle 
in a box, recovering the exact result. Our new representation of the path integral can be applied also to problems 
in which the Hamiltonian contains non-local operators (which would be the case of the Salpeter equation) and it is 
suitable both for numerical and analytical calculations, in the cases in which the limits N, N — ► oo can be calculated 
(as for a free particle in a box). 

Given the importance of path integrals in many areas of Physics, wc feel that the results contained in this paper 
could have a large number of applications. Finally, we wish to mention that it would be worth exploring the possibility 
to use the PMS in a numerical scheme to optimize convergence to the exact result for finite values of N and N. It 
also remains to explore the possibility to apply our formalism to Quantum Field Theory (QFT), which we plan to 
address in future works. 
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